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ABSTRACT 



Near-infrared reddening of extra-galactic GMCs in a face-on 

geometry 
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l/~) ' Aims. We describe the near-infrared reddening signature of giant molecular clouds (GMCs) in external galaxies. In particular, we 

examine the £j_ H and £h-k color-excesses, and the effective extinction law observed in discrete GMC regions. We also study the 
effect of the relative scale height of the GMC distribution to the color-excesses, and to the observed mass function of GMCs when the 
masses are derived using color-excess as a linear estimator of mass. 

Methods. We perform Monte Carlo radiative transfer simulations with 3D models of stellar radiation and clumpy dust distributions, 
resembling a face-on geometry. The scattered light is included in the models, and near-infrared color maps are calculated from the 
simulated data. We perform the simulations with different scale heights of GMCs and compare the color-excesses and attenuation 
5_( . of light in different geometries. We extract GMCs from the simulated color maps and compare the mass functions to the input mass 

' functions. 

Results. The effective near-infrared reddening law, i.e. the ratio £j_h/£h-k, has a value close to unity in GMC regions. The ratio 
depends significantly on the relative scale height of GMCs, f , and for f values 0. 1 . . . 0.75 we find the typical ratios of 0.6 . . . 1.1. The 
effective extinction law turns out to be very flat in GMC regions. We find the ratios of apparent ectinctions of A^/A^ = 1.35 ... 1.55 
and Aj/Ay = 1.15. The effect of the scattered flux on the effective reddening law, as well as on the effective extinction law, is 
significant. Regarding the GMC mass function, we find no correlation between the input and observed slopes of the mass functions. 
Rather, the observed slope reflects the parameter £ and the dynamical range of the mass function. As the observed slope depends on 
the geometric parameters which are not known, it is not possible to constrain the slope of the mass function using this technique. We 
estimate that only a fraction of 10 . . . 20% of the total mass of GMCs is recovered, if the observed color-excess values are transformed 
to masses using the Galactic reddening law. In the case of individual clouds the fraction can vary between ~ . . . 50%. 

Key words. Radiative transfer - Scattering - dust, extinction - ISM: clouds - Galaxies: ISM 
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• • 1. Introduction cally rules out the possibility to study them in detail using single- 

^ ' , . , . , _ , . dish observations at frequencies of common tracers of molecu- 

The dust content of galaxies has an enormous effect on their ob- ^ such afj CQ Interferometric observations can achieve the 

served properties in both local and high-redshift universe. The required resolution> but mapping large m is overw helmingl y 
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g . dust is responsible for the thermal emission at infrared wave- time consumi and sensitivity is rather pool , 

lengths, it attenuates the stellar radiation from near-infrared ^. . . . , . - . . 
. .. , j ui i • .i i - . • . ,. • , . Direct imaging of galaxies from optical to near-infared ot- 
to optical, and blocks efficiently radiation at ultraviolet wave- . ... b , . . , ,. rm .r> • -i 
i *u to, i • . - i ... c . , rp . c , . ■ fers a high-resolution view to the disttibution of GMCs inside 
lengths. The realistic description of the effects of dust is thus , . , , , , . 

. v \i c .u j them. At these wavelengths an arcsecond-scale resolution is eas- 

crucial for numerous applications ranging from the description ., , . , , , . . , . ^ ,. . , 

c , , .. . .f .. .• . ca i . ily achieved with standard imaging observations. Qualitatively, 

of galaxy evolution to the estimation of the amount of dark mat- / b r * , 

. aii. -. ■ . . .u .• .• e .u • . e a . • • the presence and distribution of GMCs can be established by 
ter. Albeit important, the estimation of the impact of dust is, in , F , _ , . . ,„,,-, 

. c .i. ii u . . • • i . i rri. - • u .u c color maps made from observations in two broadband filters, e.g. 

most of the cases, all but a trivial task. This is because the ei- . , , , b 

c . , , , . , , , .. , B — V or J — K. This approach has been used, often accompa- 

tects caused by dust depend drastically, not only on the amount . , . , . . , . , , 

c .. , . , tu . .■ i a- . u .• i ,i u • i med with image processing techniques, such as unsharp mask- 

of it, but also on the exact spatial distribution and the physical . ,. 6 , ^ , r , , , F _ 

f . ing, to disentangle the GMCs from the color maps (e.g. Regan 

"TrhapTthemost obvious manifestation of dust in galaxies ^owk & Savage HHZi™™ TrewheUaHHS 

. iu lT n , e a i i a • .u Elmegreen 1998, Thompson et al. 2004). The studies are tvpi- 

outside the Milky Way are the lanes of dark clouds seen in the „ 6 . ' r 1 . . F , , : — , ... JV 

, • i i • t f i .i_ a' ^.-u .■ c .1. cally performed in optical wavelenghts due to the higher con- 

nearby spiral galaxies. In face-on spirals the distribution of these J r , , , ,- , ■, , , , 

, . ■ . , , , , /r ^T.^ s f n v. .• i .u trast towards the dust features, and until the last decade, due to 

dusty giant molecular clouds (GMCs) follows qualitatively the , , , , ' ' 

. ,, ■ , -i a ■ i .u r^-Knr^ the lack of large near-infrared array cameras, 

stellar spiral arms. In edge-on spirals the GMCs are seen as an & 1 

optically thick layer, concentrated heavily on the midplane of The interpretation of the observed color map data is severely 

the galaxy. At the distance of nearby galaxies the sizes of GMCs hampered by the unknown relative geometry of dust and stars in 

correspond to angular diameters of a few arcseconds. This basi- the 8 alax y- The flux eimtted b y stars between a dust feature and 

the observer, as well as the star light that enters a cloud from all 

Send offprint requests to: J. Kainulainen directions and scatters towards the observer, both reduce the red- 
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dening signature of dust clouds. Thus, the description of redden- 
ing effects requires radiative transfer calculations. As the basic 
problem regarding the radiative transfer in dusty medium is im- 
portant for several fields, the studies of dusty systems have been 
made in the past by numerous authors and from several points of 
view. In the context of embedded, galaxy-like dust distributions, 
a lot of emphasis has been placed on describing the obscuration 
by inhomogenous dust which has a low volume filling factor, 
and in which large density contrasts are present (e.g. Witt et al. 
[19921 [19961 120001 Pierini et al. |2004l l. This set-up resembles the 
geometry of mixed GMCs and stellar distribution, with low col- 
umn density lines of sight dominating the light escaping from 
the system (i.e. the observed flux). It is particularly tailored to 
interpret observations of central regions of starburst galaxies, or 
high redshift star-forming galaxies. It is a common result of these 
studies that color-excesses can be highly nonlinear functions of 
optical depth and they saturate at relatively low values. The small 
scale structure, i.e. the dumpiness, significantly lowers the ob- 
served reddening compared to that caused by a homogenous dust 
distribution, and the effect of scattering is not neglible at most of 
the wavelenghts. 

In this paper we examine the obscuration effects observed 
in the regions of particularly high column density, i.e. in the 
GMCs. The configuration chosen in our simulations differs from 
the studies quoted above, as in our models the dust structures 
experience strong illumination from an extended distribution of 
stellar light surrounding the dust. The dust clouds are small com- 
pared to the extent of the stellar light distribution, and therefore 
the flux observed in on-cloud pixels is likely to be dominated by 
the emission of the stars completely in front of all the dust. The 
paper is second in series where we study the near-infrared red- 
dening properties of GMCs, which are embedded inside a disk- 
like stellar radiation field. In the first paper we showed that the 
slope of the observed GMC mass function, derived from the NIR 
color-excess data using foreground screen approximation, can 
significantly differ from the true mass function (Kainulainen et 
al. 120071 Paper I hereafter). In this paper we expand the results of 
Paper I and describe the near-infrared reddening and attenuation 
occuring in the GMC regions in detail. We examine the effective 
reddening and extinction laws in our models, with emphasis on 
how the relative scale height of the GMC distribution, affects 
their values. We also examine the shape of the observed mass 
function as a function of the parameter £ to find out, if the slope 
of the mass function can be constrained when some geometric 
parameters are known. In $2] the radiative transfer method and 
cloud models are presented. The results are presented and dis- 
cussed in |3] an d $4] In $5]we give our conclusions. 



2. Modeling 

The radiative transfer method, the basic geometry of the mod- 
els, and the cloud extraction has been described in Paper I. In 
short, the models consist of an exponential distribution of dis- 
crete clouds, which are embedded in an exponential distribution 
of stellar light. The left panel of Fig. Q]shows an example of in- 
put density distribution projected to a plane. The distributions 
of clouds and stars are characterized by their scale heights, z c 
and Zs, respectively. The ratio of the scale heights, ^ = z c /z s , is 
the main parameter describing the model geometry. The mod- 
els are viewed from a face-on angle, i.e. the inclination angle of 
zero. The dust properties in the simulations are based those of 



Draine (2003), and we use the tabulated values for Milky Way 
and Ry = 3.1, which are available on the wetQ. 

We generated 36 models, which result from choosing four 
different values for three values for the slope of the input mass 
function, a^, and three values for the low-mass cut-off of the 
mass function, M\ ow . The models and corresponding input pa- 
rameters are listed in table Q] In addition to the three parameters 
above, table Q] gives the volume filling factor, and the surface 
density filling factor in each model. The filling factor is defined 
as the number of cells with non-zero density divided by the total 
number of cells. The surface density filling factor is defined in a 
similar way, but it is calculated after projecting the model on a 
two-dimensional plane. 

The radiative transfer calculations are performed in JHK 
bands, and the resulting fluxes are used to construct simulated 
surface brightness maps. The surface brightness maps are further 
combined to color maps, namely H-K and J-H. The right panel 
°f Fig.Q]shows an example of resulting H-K color map of one 
realization of model #14. The color maps are further combined 
to Ay extinction maps using the NICER method (Lombardi & 
Alves 2001), applied to the surface brightess data. In the NICER 
method, the estimate of Ay results from the comparison of the 
colors of on-cloud pixels to the colors of nearby off-cloud pixels. 
In our models the colors of off-cloud regions are constant, and 
thus we selected one off-cloud region per model from which the 
off-cloud colors were calculated. In the selection we used the 
empty regions of J - H color maps. The individual clouds are 
identified from the extinction maps using the clumpfind routine 
'clfind2d' (Williams et al. |19941 >. The masses of identified clouds 
are calculated by summing up the extinction of all the on-cloud 
pixels, assuming the foreground screen geometry, and using the 
Milk y Way ratio ofN(H2)/Ay = 9.4 x 10 2l) cm^mag- 1 (Bohlin 
et al. |1978l >. Thus, the masses we refer to in the paper represent 
the gaseous masses of the clouds, and we express it in the units 
of solar masses, M Q . The masses are then used to construct the 
observed mass functions presented in Section[3] It is noteworthy 
that practically all the clouds that were embedded in the models 
(> 95 %) were dense enough to be detected, if they were located 
completely in front of the stellar light distribution. This means 
that there is no incompleteness in the observed mass function 
due to the family of clouds that intrinsically cannot be detected. 

The same procedure of cloud extraction was used to con- 
struct the true cloud mass functions from the inital density dis- 
tributions that were projected on a two-dimensional plane. The 
power-law model was fitted to the data to extract the slope of 
the true mass function, a true . We note that this is not exactly the 
same as a m n, which was defined to be the slope of the distribu- 
tion from which the masses of input clouds were determined. In 
other words, the true mass functions presented in Section [3] are 
not composed of the initial input masses of the clouds, but rather 
their masses after 2D projection and the clumpfind treatment. 
This was done to make sure that the cloud extraction procedure 
itself does not affect the slope of the observed mass function. 
The results of the power-law fits show that the slopes of the true 
mass functions are in good agreement with the initial ones, de- 
viating 0.2 from the value of them at most (the slopes are given 
in columns 2 and 7 of Table[T]i. 

3. Results 

The results of the radiative transfer simulations, which are the 
main results of this paper, are summarized in Table [1] The first 

1 http://www.astro.princeton.edu/~draine/dust/ 
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Fig. 1. An example of one realization of a model used in 
the simulations (model #14). Left: Column density. Right: 
Corresponding simulated H - K color map. 



five columns of TableQ]give the input parameters for each model. 
The definitions of these parameters have been given in Section[2] 
The following columns show the output properties of the models, 
atme and a b s , given in columns 6 and 7, are the exponents of the 
power-law fits to the true and observed mass functions. The 8th 
column gives the fraction of the detected clouds, i.e. the number 
of clouds extracted from the NICER extinction map divided by 
the number of input clouds. The column 9 shows the total ob- 
served mass divided by the the total input mass. The columns 10 
and 1 1 show the slope of the fit made to the frequency distribu- 
tion of color-excess pixel values (see Section U.l.ll ). 

3.1. The observables 

In the following we describe the observable characteristics of 
a typical series of simulations, and the dependency of the ob- 
servables on the model parameters. We use the models 5-8 as 
examples in most of the figures, as they represent an average set 
of input parameters. 

3.1.1. Color-excesses and the effective reddening law 

Fig. [2] shows the normalized frequency distributions of Eu-k 
and £j_h values in models 5-8. The dotted, dashed, dash-dotted, 
and solid lines are for models with £ values of 0.1, 0.27, 0.5 
and 0.75, respectively. Below -0.1 the distributions are domi- 
nated by the pixel noise. At higher color-excesses the histograms 
follow roughly an exponential distribution, with the slope obvi- 
ously depending on the value of The exponential decay of the 
color-excess values is interesting, as the initial distribution of 
cell densities, and thus the optical dephts of pixels, was a power- 
law. We fit a model logy = -Ax + B to the color-excess values 
higher than 0.2. We note that the change in the number of clouds 
per unit volume, i.e. the change in surface density, only shifts the 
distribution in vertical direction, but does not affect the parame- 
ter A of the fit. The top panel of Fig. [2] shows also the fits on top 
of the pixel histograms. The A parameters resulting from the fit 
are given in Table Q] and summarized in the lower panel of Fig. 
[2] The A parameters depend strongly on being about 10 for 
models with £ = 0.1, and shallowing with increasing values of 
In the models with £ = 0.75 the average value is ~ 4. 
The ratio of the two color-excess, E]-u/Eu-k, is especially 
interesting, as it describes the effective reddening law in the par- 
ticular geometry. We fitted a linear model to all (Eu-k, £j-h) 
pixel values of a model, and the resulting slope was ~1 in practi- 
cally all models. This is, as expected for embedded dust distribu- 
tion, much flatter than the reddening laws for Galactic dust and 
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Fig. 2. Top: Normalized histograms of Eu-k and £j_h values 
in the models 5-8 (dotted, dashed, dash-dotted, and solid lines, 
respectively). The solid lines are fits of the equation logy = 
-Ax + B to the data with E > 0.2. Bottom: Summary of the 
fitted A parameters for all the models of the study. The blue as- 
terisks correspond to Eu-k values and the red diamonds to Ej-u 
values. The error bars represent the lcr error of the fits. The hor- 
izontal lines give the error-weighted average of the data points. 
The average values are also printed into the frames. 



foreground screen geometry (in our models: Ej-u/Eu-k = 1-5, 
Draine 120031 

The ratio of color-excesses can be determined individually in 
each of the clouds. We perform this by taking the average value 
of £h~k/£j-h ratios of ten pixels with highest Eu-k to represent 
the ratio in that cloud. Fig. [3]shows the histograms of the result- 
ing ratios in models 5-8. Even though the distributions are wide, 
reflecting the large scatter in the Ej-u/Eu-k values of pixels due 
to photometric errors, they are clearly distinct. The peak of the 
distribution moves towards higher values for higher values of 
which is expected as the fraction of foreground light gets smaller 
on average. This behaviour remains hidden in the average rela- 
tion, where the distribution is dominated by the pixels with rela- 
tively low optical depths. However, it becomes detectable when 
determined separately for individual clouds. 



3.1.2. The cloud mass functions 

The masses of detected clouds are used to construct the observed 
mass functions, and a power-law model, described by the param- 
eter tt bs, is fitted to data. For the binning we use the equation 
M>ins = (2 x b) 1 ' 3 where n is the number of clouds in the his- 
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Table 1. The input parameters and observed properties of the 
models. 



Input parameters Observed properties 



# 




? 


log M low [M G ] 


// [%] 


P2D [%] 






'it ™ 






Aj-H 


t 


2.5 


0.10 


4.0 


23 


11 


- 


- 


< 1 




- ± - 


- ± - 


2 


2.5 


0.27 


4.0 


12 


11 


2.6 


2.6 


2.7 


7.3 


6.7 ± 0.75 


9.0 ± 0.33 


3 


2.5 


0.50 


4.0 


7.2 


11 


2.6 


2.2 


6.2 


11 


6.3 ± 0.78 


5.8 ± 0.30 


4 


2.5 


0.75 


4.0 


4.9 


11 


2.4 


2.2 


8.6 


22 


3.7 ± 0.49 


5.2 ± 0.30 


5 


2.0 


0.10 


4.0 


13 


14 


1.9 


2.7 


7.1 


5.8 


10 ± 0.40 


9.0 ± 0.26 


6 


2.0 


0.27 


4.0 


12 


13 


2.0 


2.7 


7.9 


7.2 


7.1 +0.26 


7.4 ± 0.24 


7 


2.0 


0.50 


4.0 


7.2 


14 


1.9 


2.0 


11 


8.8 


5.7 + 0.21 


5.2 ± 0.21 


8 


2.0 


0.75 


4.0 


5.4 


14 


2.0 


2.0 


14 


9.8 


3.7 + 0.13 


4.5 ±0.13 


9 


1.5 


0.10 


4.0 


28 


17 


1.5 


2.8 


33 


7.0 


8.4 ± 0.23 


8.3 ± 0.26 


10 


1.5 


0.27 


4.0 


14 


18 


1.5 


2.4 


34 


6.3 


7.6 + 0.25 


7.4 ± 0.11 


11 


1.5 


0.50 


4.0 


9.0 


18 


1.6 


2.0 


35 


8.3 


5.1 +0.13 


5.4 ± 0.13 


12 


1.5 


0.75 


4.0 


6.2 


17 


1.5 


2.1 


33 


11 


4.1 ±0.11 


4.5 ± 0.08 


13 


2.5 


0.10 


4.5 


28 


24 


2.6 


3.1 


8.0 


7.0 


9.5+0.85 


10 ±0.41 


14 


2.5 


0.27 


4.5 


18 


24 


2.6 


2.7 


14 


8.4 


7.0 + 0.68 


7.7 ± 0.65 


15 


2.5 


0.50 


4.5 


12 


24 


2.7 


2.8 


24 


10 


5.1 ±0.43 


6.2 ± 0.34 


26 


2.5 


0.75 


4.5 


8 


24 


2.7 


2.4 


31 


14 


3.6 ± 0.34 


4.6 ± 0.33 


17 


2.0 


0.10 


4.5 


55 


31 


2.1 


2.4 


19 


8.2 


9.2 ± 0.35 


9.5 ± 0.30 


18 


2.0 


0.27 


4.5 


38 


31 


2.0 


2.2 


21 


9.2 


6.8 ± 0.23 


7.2 ± 0.20 


19 


2.0 


0.50 


4.5 


27 


31 


2.0 


2.0 


27 


13 


4.9 ± 0.20 


5.2 ± 0.20 


20 


2.0 


0.75 


4.5 


20 


31 


1.9 


2.0 


34 


13 


4.1 ±0.13 


4.0 ±0.12 


21 


1.5 


0.10 


4.5 


27 


17 


1.5 


2.4 


17 


7.2 


9.4 ± 0.27 


8.2 ±0.19 


22 


1.5 


0.27 


4.5 


14 


15 


1.5 


2.3 


9.8 


8.5 


6.9 ± 0.25 


6.8 ± 0.20 


23 


1.5 


0.50 


4.5 


10 


16 


1.5 


2.4 


13 


8.0 


4.6 ±0.15 


4.4 ± 0.20 


24 


1.5 


0.75 


4.5 


8.0 


19 


1.5 


2.5 


19 


9.1 


3.9 ± 0.08 


3.8 ± 0.09 


25 


2.5 


0.10 


5.0 


36 


32 


2.7 


2.7 


53 


5.1 


12 ± 0.57 


11 ± 0.19 


26 


2.5 


0.27 


5.0 


20 


32 


2.5 


2.6 


56 


6.7 


8.0 ± 0.48 


9.3 ± 0.22 


27 


2.5 


0.50 


5.0 


10 


32 


2.5 


2.8 


58 


8.1 


7.1 ±0.23 


6.4 ± 0.17 


28 


2.5 


0.75 


5.0 


7.8 


32 


2.5 


2.6 


59 


11 


4.3 ±0.17 


4.5 ± 0.14 


29 


2.0 


0.10 


5.0 


37 


30 


2.0 


2.9 


68 


7.2 


11 ±0.35 


11 ±0.28 


30 


2.0 


0.27 


5.0 


14 


24 


2.0 


2.7 


91 


9.0 


7.6 ± 0.37 


7.5 ±0.31 


31 


2.0 


0.50 


5.0 


11 


30 


2.0 


2.3 


83 


10 


5.6 ±0.15 


5.5 ± 0.23 


32 


2.0 


0.75 


5.0 


7.0 


30 


1.9 


2.3 


94 


14 


4.2 ±0.11 


4.3 ±0.10 


33 


1.5 


0.10 


5.0 


28 


17 


1.4 


2.4 


71 


7.0 


8.2 ± 0.27 


8.1 ± 0.29 


34 


1.5 


0.27 


5.0 


14 


18 


1.5 


2.2 


78 


7.6 


7.8 ±0.19 


7.6 ± 0.18 


35 


1.5 


0.50 


5.0 


9.0 


18 


1.6 


1.9 


85 


8.5 


5.4 ±0.12 


5.6 ±0.16 


36 


1.5 


0.75 


5.0 


6.2 


17 


1.5 


2.1 


78 


10 


4.3 ±0.12 


4.6 ± 0.09 



" Slope of the fit to the frequency distribution of color-excess pixel 
values, i.e. logy = -Ax + B. 

b With the parameter settings of this run only a very small fraction of 
the clouds were detected. 



togram. The lowest mass bin which is taken into the fit is selected 
individually for each model by eye. The slopes of both true and 
observed mass functions are listed in columns 6 and 7 of Table 
[TJ The appearance of the mass functions is similar to those de- 
scribed already in Paper I (Fig. 2). The dynamical range of the 
observed mass function is approximately two decades, starting 
from ~ 10 3 ' 5 M Q and terminating at 10 5 5 M . In some mod- 
els the power-law does not represent the observed mass function 
well over the whole observed dynamical scale. Rather, the mass 
function has a break-point at M ~ 10 4 5 M , below which the 
mass function flattens, or is completely flat. This typically hap- 
pens in the models where both £ and the low-mass cut-off of the 
mass function, M\ ow , are high. 

Fig. |4] shows the relation between the true and observed 
slopes in all the models. The observed slopes vary between 1.8 
and 2.9, and there is no correlation between the true and ob- 
served slopes. The models with flat input slopes are observed 
to have similar range of slopes as the models with steep input 
slopes. In some particular parameter settings, such as models 



marked with purple diamonds (models 27-28, 31-32, and 35- 
36), the observed slopes follow loosely the true slopes if the £ 
is high. In these models the M\ ovl parameter is 10 5 M . The ob- 
served slopes correlate slightly with the parameter as on av- 
erage, the models with low £ values tend to show steeper slopes 
(<*obs ~ 2.5) than the models with high £ values (a Q \, s ~ 2). As 
the correlation between a trae and ar b s is very poor, it is not pos- 
sible to estimate atme by using simply a measurement of a \, s . 

3.1 .3. The sizes of the clouds 

The sizes of the detected clouds, i.e. the number of on-cloud 
pixels, were compared to the sizes of the corresponding input 
clouds. The observed sizes are smaller than input sizes with an 
average factor of 0.85. The ratio remains remarkably constant, 
as 90 % of all clouds in all the models have ratios between 
0.7 ... 1.0. Thus, if the cloud is detected using near-infrared ex- 
cesses, it is likely to be detected close to its full extent, and only 
the most diffuse outer edges remain undetected. 
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Fig. 3. The histograms of the £j_h/£h-k ratios of all detected 
GMCs in the models 5-8. The ratio is measured at the point of 
maximum £h-k in each cloud. The approximate peak values 
corresponding the £ values of 0.1, 0.27, 0.5, and 0.75 are 0.5, 
0.6, 0.9, and 1.1. The number of the model and the correspond- 
ing £ value are marked into the frames. The dotted lines show 
the ratio if the scattered flux is neglected (see section [3"3| >. 



3.1.4. The apparent extinction 

As the absolute value of background flux in our models is known 
(the background flux is constant and the positions with zero op- 
tical depth are known), we can examine the apparent extinction 
in the models. The apparent extinction is defined as: 



F(A) 

F] 



(1) 



where is the unattenuated background flux. The values of 
apparent extinction are typically in the order of a few tenths of 
a magnitude in our models. We estimate the wavelength depen- 
dency of the apparent extinction by calculating the ratios A^/A^ 
and A^/A^ at the positions of each extracted cloud. These ratios 
describe how the extinction depends on wavelength in the partic- 
ular geometry, and thus the curve defined by them is often called 
an effective extinction law. We calculate the ratios in each cloud 
by taking the average of ten pixels where the £h-k is highest. 
The resulting histograms are shown in Fig. The distributions 
of Aj^/A^ values have their maxima between ~ 1.35 .. . 1.55, de- 
pending slightly on the value of The distributions of AJ/A^ 
have their maxima only slightly above unity, typically at ~ 1.1. 
The ratios are clearly lower than in the case of foreground screen 
extinction (in our models: tj/th = 1-58 and th/tk = 1.65, 
Draine 120031 . 



3.2. The amount of undetected mass 

As was already found in Paper I, the masses of the clouds that 
are detected can be drastically smaller than their true masses. 
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Fig. 4. The slope of the observed mass function as a function of 
the slope of the true mass function in all models of the study. The 
different panels correspond to different values of £ parameter. 
The red plus signs correspond to Mj ow = 4.0 (models 1-12), blue 
asterisks to M\ ovl = 4.5 (models 13-24), and purple diamonds to 
Miow =5.0 (models 25-36). 



As an example, Fig.|6]shows the comparison of the true and ob- 
served masses of clouds in model 8. The scatter in the relation 
is very high. Essentially, a cloud can be observed with any mass 
smaller than its true mass with a significant probability, down to 
the detection limit. In addition to the mass that is missed due to 
the underestimation of the masses of detected clouds, a fraction 
of clouds always remains completely undetected. This is mainly 
due to the population that is located on the further side of the 
galaxy. As the flux of the on-cloud pixels is completely domi- 
nated by the foreground flux, these clouds do not produce signif- 
icant reddening features to the color maps. We calculate the total 
amount of undetected mass by summing up the observed extinc- 
tion maps, and comparing them to the total sum of input column 
density maps. The resulting ratios, reported in column 10 of table 
[TJ are obviously very dependent of what is the dynamical range 
of input masses. Thus, the ratios given should only be taken as 
estimates of the order of undetected mass. 

Following the treatment in Paper I, we also calculated the 
fraction of detected clouds as a function of true mass by com- 
paring the positions of clouds detected from input column den- 
sity maps to the positions of clouds detected from the NICER 
extinction maps. The behaviour of the resulting completeness 
function was described in paper I. The larger set of models in 
this paper allows the examination of £ dependence of the com- 
pleteness function. Fig.|6]shows an example of the £ dependence 
in models 5-8. The decrease in the completeness below 10 5 5 M 
is stronger for smaller values of In the models with £ = 0.1 ba- 
sically no clouds with M < 10 4 5 M are detected. In the models 
with \ = 0.75 clouds with M = 10 4 M Q still have a probability 
of a few percent to be detected. Above 10 5 5 M the detection 
probability is slightly higher for the models with lower £. This 
is simply due to the fact that the amount of clouds for which the 
fraction of foreground flux is high is larger in models with high 
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Fig. 5. Top: The histogram of the ratio of apparent extinction, 
Aj^/A^, evaluated in each on-cloud region. The apparent extinc- 
tion is as defined by Eq.Q] The data are for models 5-8, with the 
model number and the corresponding £ values are marked into 
the frames. The dotted lines show the corresponding histograms 
if the scattered light is neglected, (see section [3~3l Bottom: The 
same for Aj/A^. 



The color of foreground light dominates the colors observed 
on these clouds over the small fraction of reddened background 
light. 

3.3. The significance of the scattered light 

During the simulations the scattered and emitted fluxes were reg- 
istered separately, so we can explore the role of the scattering in 
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Fig. 6. Top: The observed masses of clouds as a function of 
the true mass of clouds in model 8. The observed masses are 
calculated from the £h-k values by assuming the foreground 
screen geometry, and Galactic gas-to-dust ratio. The dotted line 
shows the one-to-one ratio. Bottom: The fraction of detected 
clouds, i.e. the completeness function, in models 5-8 (dotted, 
dashed, dash-dotted, and solid lines). The masses are given in 
Solar masses. 



the effective reddening and attenuation laws by using only the 
emitted fluxes. Figs. [3]and[5]show the corresponding reddening 
and attenuation laws with dashed lines. As the scattered light is 
not present the dispersion in histograms is much smaller. The 
peaks are clearly shifted towards higher ratios in both reddening 
and attenuation data. Thus, the scattered flux has a significant 
effect on both values. 

The level of the scattered flux with respect to the total ob- 
served flux is shown in Fig. [7] To better visualize the data, the 
figure shows the probability density of data points smoothed 
with a gaussian kernel. The data is shown for the model 5, and 
the overplotted lines show the median value of the relation in 
models 5 and 8. The increase of the fraction of the scattered flux 
is close to linear up to Ay ~ 3 . . . 5'". At higher optical depths the 
fraction saturates quickly, reaching its maximum at Ay ~ 10 m 
in J and H bands. In K band the fraction continues to increase 
slowly until Ay ~ 30™. The value of £ has a relatively small, 
yet detectable, effect to the median value of the relation. The 
median value decreases slightly with increasing which results 
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4. Discussion 
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Fig. 7. The density contours of the fraction of scattered light in 
model 5 at Top: /, Middle: H, and Bottom: K bands. The frac- 
tion refers to the scattered flux in the pixel divided by the total 
flux in that pixel (i.e. not by the total unattenuated flux). The 
solid lines show the median value of the relation for different Ay 
values in model 5 = 0.1). For comparison, the dashed line 
shows the median ratio for model 8 = 0.75). 



from the average radiation field experienced by a cloud being 
stronger when ^ value is low. 

In some studies an estimation of the optical depths along 
the line of sight has been attempted by using measurements of 
apparent extinction at two or more wavelengths (e.g. Howk & 
Savage 119971 12Q001 >. In the approach the scattered flux is ne- 
glected, and the resulting equations for radiative transfer become 
possible to solve. We applied this method to our models, but 
the equations failed to converge for most of the clouds. This 
failure was clearly due to the scattered flux, as the same pro- 
cedure using only emitted fluxes indeed yielded improved esti- 
mates of optical depths. We also made an effort to use the data 
presented in Fig.Qto solve numerically the actual optical depth. 
Typically, this operation multiplied cloud masses with a factor 
varying from 1 . . .5. Even though on average the cloud masses 
were thus improved, the effect was highly variable, and in every 
model the form of the mass function was completely destroyed. 
Considering that in our models the uncertainty caused by un- 
known background flux is neglible, we conclude that this cor- 
rection is not feasible to gain an improved estimate of the mass 
function. 



As has been well established in the previous studies of the sub- 
ject, the large scatter and strong saturation of colors as a function 
of the total optical depth are dominating the simulations. These 
effects, originating from the varying relative geometry and in- 
homogeinity of dust distribution, affect both the level of extinc- 
tion, and its wavelength dependency (i.e. the effective extinction 
law). The knowledge of realistic effective extinction law is cru- 
cial for numerous applications, where the observed signal from 
an object or region of interest needs to be corrected for the ef- 
fects of dust. For example, for this purpose Calzetti et al. ( 1997 
2000 2001) have derived an empirical attenuation laws for cen- 
tral regions of starburst galaxies. Likewise, theoretical attenua- 
tion curves have been calculated for several geometries of inho- 
mogenous dust and stars (e.g. Pierini et al. 2004 and references 
therein). In this paper we particularly examine the on-cloud pix- 
els rather than the integrated light over a region covered with a 
mixture of GMCs and low column density lines of sight. We em- 
phasize how the reddening properties of discrete cloud regions 
can vary, even in the fixed face-on geometry, when embedded 
GMCs are considered. 

As we are considering only the on-cloud pixels, the effects 
due to the scattered flux undoubtly increase. In the most massive 
clouds, which are most likely to be detected in extinction maps, 
the fraction of the scattered flux to the total flux of a pixel can 
go up to several tens of percents (Fig.|7J. The scattered flux has 
a significant effect to the scale of color-excesses, dampening the 
overall color-excesses level with a factor of 1 .2 ... 1 .7. 

Our simulations predict the effective reddening law, con- 
structed from the values of the on-cloud pixels of detected 
clouds, to be close to unity. This result is valid only for the se- 
lected geometry, i.e. for a face-on galaxy. Even though the same 
average relation was found in all the models, we showed that 
when each cloud is examined separately, a £ dependence can 
be found (Fig. This demonstrates how a relatively straight- 
forwardly measured observable reflects the geometry of the ob- 
served region despite the complicating effect of the scattered 
flux. In practice the determination of Ej-h/^h-k ratio for a sin- 
gle cloud is rather uncertain due to the photometric accuracy and 
generally small values of £j_h and £h-k, an d thus a sample of 
hundreds of clouds is likely to be needed to gain an estimate 
of the £ value of a face-on galaxy. It is interesting that even in 
the models with lowest average foreground emission, i.e. models 
with £ = 0.75, the effective reddening law remained very flat. 
Thus, if the overall conditions do not differ significantly from 
the simulations, high £j_h/£h-k ratios should not be observed 
unless the foreground emission in the observed region really is 
neglible. This means that if a cloud region is observed to have 
a high £j_h/£h-k ratio, it can be argued that the foreground 
screen approximation is valid for it. 

In the case of an embedded cloud the effective reddening law 
does not relate the observed color-excesses directly to apparent 
extinction (as it would in the case of a foreground cloud). Thus, 
the effective extinction law cannot be inferred from those data 
alone. We evaluated the effective extinction law by determining 
the apparent extinction of on-cloud pixels, which is trivial to do 
in our models. We do note that this is not typically the case with 
spiral galaxies, as the knowledge of absolute background flux 
to a good precision would require a detailed model of the dis- 
tribution of stellar light. The typical ratios of AJ/A^ = 1.15 and 
A^/AjJ, = 1.35 ... 1.55 were found (Fig. The ratios are similar 
to the values derived in a pair of occulting galaxies by Berlind 
et al. 119971 We find this similarity quite convincing, as the re- 



8 



J. Kainulainen et al.: Near-infrared reddening of extra-galactic GMCs in a face-on geometry 



gions sampled in Berlind et al. are relatively small, not much 
larger than the sizes of individual GMCs. Thus, the study might 
represent the conditions described by our simulations relatively 
well. 

In addition to the description of reddening effects, we exam- 
ine the mass functions derived from the color-excess data, and 
thus expand the work done in Paper I by increasing the number 
of models. The main emphasis regarding the mass functions is to 
examine the observed slope at different values of In the context 
of cloud mass functions, the direct consequence of the satura- 
tion of color-excesses is that the observed mass function is con- 
fined to a narrow dynamical range compared to the original mass 
function. A comparison of models which are otherwise similar, 
but have different £ values, shows that the confinement depends 
slightly on The difference in the high-mass end between the 
models with £ = 0.1 and £ = 0.75 is typically logM ~ 0.4 M , 
while the low-mass end is at log M ~ 3.5 M in all the models. 

In most of our models the observed slope of the mass func- 
tion differs at least 0.3 from the slope of the true mass function. 
As is shown if Fig. |H there is practically no correlation between 
the true and observed slopes. The observed slope depends ob- 
viously on the parameter as in the models with ^ = 0.1 the 
slopes are generally > 2.5, and in the models with £ = 0.75 they 
are closer to 2. However, £ clearly is not the only parameter af- 
fecting the slope. Whereas the models with M\ ovl — 10 4 M Q show 
practically constant observed slopes for all true slopes (blue plus 
signs in PigSJ , the slopes in the models with M\ ow = 10 5 M do 
increase with the input slope in the cases where £ is large (green 
asterisks). The reason for this is most likely that the dispersion 
between the true and observed masses of low mass clouds is 
higher than in high mass clouds. Further, when the £ value is 
high the reddening signatures are stronger and the cloud defini- 
tion better. When both £ and M\ ow are high, the relation between 
the true and observed masses is significantly better than for low 
£ and Mi ow . Thus, if there is no reason to assume a truncation on 
the low-mass side of the mass function, the presence of lowest 
mass clouds makes the slope of the observed mass function quite 
insensitive to the true slope. However, if the truncation occurs at 
significantly lower masses than 10 5 , only a relatively small num- 
ber of detected clouds will be affected by it, and thus the cut-off 
should be meaningless regarding the observed mass function. 

The near-infrared color-excess values, used in conjunction 
with the simple foreground screen approximation, do not pro- 
vide a feasible method to determine the mass function of extra- 
galactic GMCs. However, the information carried by color- 
excess data regarding the geometric parameters remains interest- 
ing as such, and its sensitivity to dust features makes it a power- 
ful tracer of GMCs despite the poor correspondance to real mass. 
Also, the results of this study are restricted to a face-on geome- 
try where the most important factor, the fraction of foreground 
flux, is high and more importantly unknown. Several situations 
can be imagined, where the distribution of GMCs with respect 
to the main illumination source is much better constrained (e.g. 
nuclear-centered rings, dust structures in highly inclined sys- 
tems, etc.). The use of NIR color-excess as a mass tracer in these 
cases may prove to be feasible, but that remains to be confirmed 
with simulations specific for the geometry in question. 

5. Conclusions 

We have examined the near-infrared reddening properties of 
extra-galactic GMCs with Monte Carlo radiative transfer mod- 
els. Our models resemble the face-on geometry and correspond 
to the dimensions of the most nearby spiral galaxies. We have 



investigated the effect of the relative scale height of GMCs, 
to the observed color-excesses and the mass spectrum of GMCs. 
We also derive the values of effective reddening law, £j_h/£h-k, 
and the effective extinction law for our models. The main con- 
clusions of our simulations are as follows: 

1. The effective NIR reddening-law of GMCs in a face-on 
geometry is predicted to have an average slope close to 
unity. The observed £j_h/£h-k ratios of on-cloud regions 
depend heavily on the exact geometry of the system, which 
in the context of a face-on galaxy is reflected by the value 
of From the simulations we find the average ratios of 
£ , j_h/£'h-k=0.6, 0.7, 0.9, and 1.1 for the £ values of 0.1, 
0.27, 0.5, and 0.75, respectively. The ratios towards spatially 
separate GMCs should be distinguishable within the typical 
photometric errors, and thus the value can be used as an in- 
dicator of line-of-sight geometry. 

2. The apparent extinction of GMC regions shows very flat 
wavelength dependency in /, H, and K bands. We find that 
typically AJ/A^ -1.15 and Ay ~ 1.35 . . . 1.55, the latter 
depending slightly on 

3. The comparisons of the models with and without scattering 
suggests that the scattered flux has a non-neglible effect to 
the color-excesses and the apparent extinction of GMC re- 
gions. 

4. We extend the results of Paper I, confirming that near- 
infrared color-excess is an unrelieable measure of the col- 
umn density of GMCs. If £h-k values are transformed to 
masses assuming the Galactic reddening law, on average a 
fraction in the order of 10-20% of the total mass is recovered. 
In the case of individual GMCs the percentage can vary alot, 
being basically anything between . . . ~ 50%. 

5. The mass function constructed by using the color-excesses 
and the Galactic reddening law is sensitive to the parame- 
ter The models with smaller £ values tend to have steeper 
(> 2.5) slopes than models with high £ values (> 2.0). The 
dynamical range of the underlying mass function also mod- 
ifies the observed slope. There is a strong degeneracy in the 
observed slope with respect to the parameters a true , and 
Mi ow . As none of these parameters is observationally well 
constrained, it is not possible to determine a true based only 
on a measurement of a b s . 
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